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THE  EFFECTS  OF  ELECTRON-NEUTRAL  COLLISIONS 
ON  THE  INTENSITY  OF  PLASMA  LINES 


I.  INTRODUCTION 


The  ionosphere  is  a  region  of  weakly  ionized  plasma  essentially  in  thermal  equilibrium.  When 
incoherent  radar  is  used  to  probe  the  ionosphere,  VHF  radiation  is  scattered  off  density  inhomo¬ 
geneities  at  various  altitudes  and  the  spectral  character  of  the  scattered  wave  is  examined.  If  the  inho¬ 
mogeneities  have  wavelengths  longer  than  the  Debye  length,  the  scattered  spectra  display  sharp  reso¬ 
nances  near  the  plasma  frequency.1  The  intensity  of  these  resonances,  called  plasma  lines,  depends  on 
the  electron  density  distribution  and  is  enhanced  by  the  presence  of  nonthermal  electrons.2  These 
nonthermal  electrons  may  be  photo  electrons  excited  by  solar  electromagnetic  radiation  or  they  may  be 
secondary  electrons  excited  by  energetic  particles  such  as  electrons,  protons  or  heavier  ions  deposited  in 
the  ionosphere.  The  shape  and  intensity  of  the  plasma  lines  are  further  altered  by  Landau  damping  and 
collisional  broadening. 

The  intensity  of  a  plasma  line  has  been  calculated  by  Yngvesson  and  Perkins’,  who  included  the 
effect  of  Coulomb  collisions  between  electrons  and  ions  as  well  as  the  effect  of  Landau  damping.  Their 
result  may  be  expressed  in  terms  of  an  effective  temperature  for  the  plasma  line. 


fm  +  fn  +  Xa 

fm  +  Lp  + 


The  terms.  fn  and  Lp ,  represent  excitation  and  damping  of  plasma  waves  due  to  nonthermal  electrons. 
The  term,  describes  Landau  damping  while  describes  effects  of  electron-ion  collisions.  The 
term  /,„  incorporates  an  assumption  that  the  distribution  of  background  thermal  electrons  is  basically 
Maxwellian.  T„  of  Eq.  (1)  may  be  related  directly  to  the  power  of  the  returned  radar  signal. 


The  measurement  of  plasma  lines  is  of  current  interest  due  to  both  the  availability  of  accurate 
observational  data4  5  and  the  development  of  sophisticated  theoretical  and  computational  techniques''  "1 
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for  calculating  the  ionospheric  electron  distribution  function.  For  example,  comparisons  of  model  pred¬ 
ictions  with  recent  Chatanika  data'  have  shown  good  agreement  at  altitudes  high  enough  to  allow  the 
neglect  of  electron-neutral  collisions. 

The  purpose  of  the  present  work  is  to  extend  the  theoretical  predictions  of  Yngvesson  and  Per¬ 
kins3  to  include  self-consistently  the  effects  of  electron  collisions  with  neutral  species  which  must  be 
included  at  lower  altitudes.  The  approach  we  have  taken  is  to  apply  Fluctuation  Dissapation  theory  ,0. 
representing  collisions  by  a  BGK  collision  term  in  the  Boltzmann  equation.  Since  Fokker-Planck  theory 
is  only  appropriate  for  Coulomb  collisions,  we  cannot  use  the  same  approach  as  Perkins  and  Salpeter. 

In  (he  next  section  we  discuss  the  representation  of  collisions  in  the  Boltzmann  equation.  The 
following  sections  outline  the  calculation  of  the  spectral  density  function,  which  is  integrated  across  the 
plasma  resonance  to  yield  the  intensity  of  a  plasma  line.  The  result  is  shown  to  compare  well  with 
effective  temperatures  derived  from  Chatanika  data. 

II.  THEORETICAL  REPRESENTATION  OF  COLLISIONS 

Perkins  and  Salpeter  applied  a  Fokker-Planck  collison  term  in  order  to  be  able  to  model  small- 
angle  deflections  for  electron-ion  collisions.  They  derived  a  collision  frequency  for  these  long-range 
Coulomb  collisions,  which  had  the  same  form  as  a  typical  momentum  transfer  collision  frequency.  This 
frequency  contains  kinetic  information  about  local  collisional  processes,  but  also  contains  an  implicit 
assumption  that  the  distribution  function  for  electrons  is  basically  Maxwellian.  Since  electron-ion  colli¬ 
sions  in  the  ionosphere  principally  involve  particles  of  average  energy,  these  collisions  cannot  produce 
an  enhancement  of  the  electron  density  fluctuations,  but  only  damping. 

In  order  to  understand  what  kind  of  approximations  can  be  made  in  representing  different  kinds 
of  collisions,  it  is  important  to  know  the  length  scales  and  timescales  which  dominate  physical  processes 
in  various  regions  of  the  ionosphere.  The  collision  time,  r  -  I/p.  represents  the  amount  of  time 
needed  for  an  electron  to  experience  a  significant  deflection.  That  is.  such  a  deflection  results  from 
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many  small  angle  collisons  in  the  case  of  electron-ion  collisions,  whereas  a  single  electron-neutral  colli¬ 
sion  would  produce  on  average  the  same  deflection.  Both  electron- neutral  and  electron-ion  collision 
times  are  much  larger  than  the  reciprocal  electron  plasma  frequency,  <u“'. 

<<*',.„<<  ti>p  :  <  150  km 

ve„  <<  Vei  <<  <op  :  >  200  km. 


The  radar  only  detects  waves  of  a  particular  wavelength,  typically  23cm  for  recent  Chatanika 
experiments.  This  is  long  compared  to  typical  electron  Debye  lengths  which  are  of  the  order  of  a  few 


centimeters.  Typical  length  scales  compare  with  the  collisional  mean  free  paths  / 


ykT 

m 


1/2 


1 

— ,  as 

V 


A*  <  /,»  <  <  A.  100  km  <  r  <  110  km 

Ao,.  <  A  <  <  /„,  r  >  120  km. 

This  relationship  suggests  that  one  could  even  use  fluid  theory  to  describe  electron-neutral  collisions  in 
the  E-region  since  the  E-region  plasma  is  not  sensitive  to  details  of  individual  collisions.  Kinetic  theory 
is  required  for  the  F-region  where  collisional  mean  free  paths  are  greater  than  the  radar  wavelength. 
Hence,  we  have  chosen  a  BGK  representation  of  the  collision  term  in  the  Boltzmann  equation.  The 
BGK  collision  term  is  particle  conservative  and  is  designed  for  electron-neutral  or  ion-neutral  colli¬ 
sions."  It  requires  that  a  collision  frequency  be  provided  which  contains  information  about  individual 

collisions  and  represents  the  frequency  at  which  significant  deflections  occur  in  the  direction  of  particle 
motion. 


Note  that  the  Perkins  and  Salpeter  electron-ion  collision  frequency  provides  an  adequate  represen¬ 
tation  which  can  sensibly  be  applied  in  this  model.  Recall,  their  value  does  not  indicate  the  frequency  of 
small-deflection  collisions  but  rather  an  averaged  value  based  on  a  collection  of  small  angle  collisions 
sufficient  to  produce  a  significant  deflection.  Since  we  will  be  assuming  that  the  electron  distribution 
function  is  basically  Maxwellian  (as  was  assumed  in  the  derivation  of  Eq.  ( 1 )  )2,  we  can  represent  the 
effect  of  electron-ion  collisions  by  a  BGK  term.  This  treatment  is  not  typical  but  may  be  used  for  cal¬ 
culating  the  intensity  of  a  plasma  line  for  which  collisions  are  a  minor  effect.  We  show  later  that  our 
results  reduce  to  the  Fokker-Planck  results  when  electron- neutral  collisons  are  neglcgible. 
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The  BGK  representation  of  the  Boltzmann  collision  term  has  the  form 

A'i(x.r) 


iL 

3/ 

*  V 

c 

/i(x.  v.  /)  - 


N  „ 


/o<v) 


(2) 


where  /0  is  the  ambient  election  distribution  in  velocity  space,  J\  is  a  perturbation  on  this  distribution, 
and  .V0  and  .Vt  are  integrals  of  Jo  and  J\  respectively,  taken  over  velocity  space.  A  Boltzmann  collision 
term  must  be  provided  to  represent  each  kind  of  collision  of  electrons  with  other  species.  The 
appropriate  collision  frequencies  have  been  derived  by  Banks13  and  Itikawa14  for  electron-neutral  colli¬ 
sions.  The  collision  frequencies  do  not  depend  on  the  electron  density  and  are  so  slowly  varying  in 
space  and  time  that  they  can  be  assumed  to  be  constant  except  for  a  weak  altitude  dependence.  [The 
quantity  has  a  value  of  10  km  at  an  altitude  of  105  km  for  the  dominant  electron-neutral 

collisions.  This  value  is  much  larger  than  the  wavelengths  of  interest  to  us.] 


A  total  collision  term  described  by  Eq.  (2)  can  hence  be  written  using  a  total  collision  frequency 

V  =  X  +  v»i 

/ 

where  the  summation  is  taken  over  all  relevant  electron-neutral  collisional  interactions. 

III.  CALCULATION  OF  THE  SPECTRAL  DENSITY.  «(*,«) 


Given  a  suitable  total  collision  frequency  vh  the  first  order  Boltzmann  equation  for  species  /  is 
dFu{r.  v.  r)  9Fu(t,  v.  r)  Q  ^  3F0/(v> 


3/ 


+  v 


dr 


+  ■  — ~ - -  -  rlF^lr.  v./>  -  h,F()/(v)] 


m  3v 


(4) 


with 


flt/tr.f)  —  J*  </vf|/lr.  v./). 


(5) 


We  assume  Fni(v)  is  a  zeroth  order  distribution  function  normalized  so  that  the  integral  of 
taken  over  velocity  space  is  unity.  F\t  (r,  v,r)  is  a  first-order  spatial  perturbation  normalized  by  the 
same  constant.  Similarly,  we  define  a  perturbation  electric  field  E|  =  -Vd.  where  6  is  the  electrostatic 
potential.  so  that 

V  •  E,  -  £  J  ihFu  (6) 
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The  summation  is  taken  over  species  /,  representing  electrons  and  ions,  of  charge  q.  We  then  have 

—  -4ir(Zen\,  —  enu.)  (7) 

In  order  to  determine  the  spectral  densities,  /?|/(k.cu),  we  calculate  the  Fourier-Laplace 
transforms  of  (7)  and  (4): 

—k2(t>  —  —  47r[Ze«|,(k.a»)  -  e«i,.<k.fc»)I  (8) 


-Fi,(k.v.O)  +  fiat  +  —  /'(k  •  v)]Fu(k,  v,<u) 


-  [Zenu(k.m)  -  e'/i u- tk. a> )1  k  • 

mk-  dv 


+  t>inu(k,o>)F0i  -  iiuik.O)  F0l 


<9) 


The  first  term  in  Eq.  (9)  appears  as  a  constant  of  integration  because  we  are  performing  a  Laplace 
transform  rather  than  a  Fourier  transform  in  time.  When  we  divide  Eq.  (9)  by  (w  -  k  ■  v  -  it',)  and 
integrate  over  v,  the  first  term  is  expressed  in  terms  of 


«,,.(k, 0)  -  £  e'k  r  WS(v-  v(0)). 

/-i 


and  (9)  becomes12: 


«,,( k.0) 


k 


BF0„ 


io  —  k  v  —  ivt. 


+  />»!,. (k.<u)  -  +  [Zgg[,(k.w)]-eff[,.(k.<«»)l  f  - - — 

mk-  J  <•>  —  k  •  v  —  i 


k  ■ v  —  iv 


F{)l. 


kv  — n\ 


We  define: 


r°°  dvFoi  . 

£Mk.w)  =  iv,  j  - r - — •  and 

a.  —  k  •  v—  n>i 


kv-  iv, 


C'<kw)  “  TTT^T J—  lh  m,kl  (u-k  y-iv. 


(10) 


(11) 


(12) 


(13) 
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Then  Eq.  UP  for  electrons  and  an  analogous  equation  for  ions  may  be  combined  to  eliminate  n,M 
yielding: 


Hi,,(k.w) 


-  / 


1- 


C 


€  (k.cu) 


1  ik  r,(0) 

•  A.  II 

1  +  a  hl 

1 

1 

1 

© 

>• 

1 

3 

zc. 

1 

•X  ik  r  <01 

1 

_  D_ 

€  (k.a») 

(1  +  D,) 

2.  e 

/-i 

cj  —  k  •  v, (0)  —  if, 

iv, 

with 


(14) 


«(k.to)  —  1  +  C,,(k,tu)  +  C,  (k,  a> )  (15) 

This  value  will  be  used  in  the  next  section  to  calculate  the  effective  temperature  associated  with  a 
plasma  line. 

IV.  CALCULATION  OF  THE  EFFECTIVE  TEMPERATURE.  T. 


As  in  the  treatment  by  Yngvesson  and  Perkins3,  we  define  a  spectral  density  function  S(k.w) 
which  describes  the  spectral  amplitude  as  a  function  of  wavenumber  and  frequency. 


2v,  <|//,,.(k.w)|:> 

S (k.w)  -  lim 


i  —  oc  1 


H..0 


(16) 


The  intensity  associated  with  a  resonant  peak  is  proportional  to  the  effective  temperature.  T„  of  Eq.  (1), 
and  is  simply  5(k.o»)  integrated  over  the  frequency  range  about  the  plasma  resonance. 

When  the  value  of  «i,,(k.w)  is  substituted  from  Eq.  (13).  the  function  S(k.u>)  becomes 

Slk.to)  -  2 1 < I  +  C,)/e|‘  B,.  +  2ZlC,./e|:fl,.  (17) 


where 


B, 


ll  +  D,\2 


</v/0,(v) 


\D,\2 


(w  —  k  •  v)’  +  vr 


(18) 


The  distribution  function  /,„(v)  describes  species  /  in  the  thermal  plasma  and  it  is  best 
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Choosing  the  Maxwellian  distributions: 

./;,0(v)  —  exp  (-v2/a2)/(ira2)il2,  a  =  (2kTJ  me)l/2, 

—  exp(—  v2/  b2)/(irb2)}/2,  b  —  (2k  TJm,)1'2, 

the  previous  coefficients  become, 

/^(k.to)  -  [2exp(— y2)  jro  exp (p2)dp  +  iir[ 12  exp(— >v2)] 

Djik.oi)  -  ~  [2exp(->',2)  J0  exp(p2)dp  +  in'12  exp(-j,2)]. 

where 


ye  =  (w  -  in,.)/ ka.  y,  =  (m  -  iv,)/kb. 


Bc(k.ai) 


1  r +°° 

ni/2ka\\  +  DeY  J-» 


,  v_  exp[-(v/fl)2] 
1  a  (Vv,.  -  v/ a) 


lD,l2 

r,.|l  +  D,  I2 


‘k  |i  +  ^  i2  Im|2exp(-v>(,2)  jn  exp (p2)dp  +  /irl/2exp(-y,2) 


\De\2 


v,\\  +  A.  I2 


fl,(k,o>) - — — i - 

w1  2M|f  +  D, |2  {•'-•»  [  b 


exp{-(v/6)2l 


(.v,  -  v lb) 


ID  I2 


11  +  ol- 


Cfk.w)  =  —  [  1  -  2v,,  exp(— v,,2)  fg  '  exp  (p2)  dp  -  inl2y,  exp(-r,:)j 


C,( k.a»)  -  | y  “  —  j  [  1  -  2p,  exp(- v,2)  JQ  exp(p:)dp  -  it rl,:v,  exp<-»,2)]. 


where  a  —  (AA.0)-1  and  AD  is  the  Debye  length.  Incorporating  these  values,  Eq.  (17)  reduces  to 


S,.(k, «»,«',.)  - 

ka 


Im[fir,  2exp(— vr2)  +  2exp(-y,.2)  exp(pr) —  —  \D,  |2 


1 1  +  De  4-  a2(Rw{y,.)  -  i7w(_v,.))|; 


Since  w  a  »  ka,  that  is  a2  >  1,  we  define  w  by 
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0)r 

k 

^  k*<arK  7> 

irmeut r 

,  ,  o(v) 

(  2 u>-vk2 

^  8v 

V  -  -r-  trto/ 

(34) 


Since  nonthermal  electrons  constitute  only  a  minor  proportion  of  available  electrons,  we  represent  them 
as  a  perturbation,  fp(y).  on  the  Maxwellian  distribution  of  thermal  electrons,  That  is. 


y;,0<v)  Jm(\)  +4<v).  (35) 

We  describe  the  intensity  in  terms  of  an  effective  temperature.  Tp ,  analogous  to  the  value,  T  defined 
by  Perkins  and  Salpeter1, 


inv/kT  \  =  —  d  \n  feo(\)/ dv 


(36) 


Note  Eq.  (34)  reduces  to  /„ 


2a -T(, 


when  v  -  0  and  T • ,  reduces  to  Te  for  ./;,0  —  Jm .  We  define  an 


effective  temperature,  71,  such  that 


lr  -  ~r=-  for  ^0. 

2a' I,, 


(37) 


Also  we  define 


_  *jk7; 

X  = - J  v 

Incorporating  Eqs.  (35),  (37)  and  (38)  into  Eq.  (34)  we  may  write 


(38) 


2 

'  V 


Tn  ./;„(«,/«)  +  J'JuiJk)  +  x 


./;„(«,./*)  +  Lp(u>r/x)  +x 


(39) 


where 


Lr  = 


—  —  3/,<v)/dv 
<itr  me 


(40) 


The  term  Lr  represents  a  loss  of  energy  due  to  changes  in  the  nonthermal  electron  distribution. 


The  effective  collison  frequency,  v,  given  in  Eq.  3,  is  the  sum  of  electron-neutral  and  ion-neutral 
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collision  frequencies.  Our  calculation  of  the  effective  electron  temperature  associated  with  a  plasma 
line.  r„  in  Eq.  (38).  reduces  to  Eq.  (1)  if  the  electron-neutral  collisions  are  negligible. 

V.  CONCLUSION 

Electron-neutral  collisons  are  only  important  in  the  lower  ionosphere,  where  the  terms  included  in 
Eq.  (38)  vary  as  shown  in  Fig.  1.  Although  the  electron-neutral  collision  terms  are  small,  they  can 
have  a  significant  effect  on  the  effective  temperature  at  low  altitudes.  This  is  seen  in  Fig.  2.  The  solid 
curves  are  calculated  using  Eq.  (38),  while  the  dashed  curves  have  been  calculated  neglecting  electron- 
neutral  collisions.  Experimental  data  from  Chatanika8,  plotted  on  the  same  graphs,  reflect  a  much 
better  correspondence  to  the  calculations  which  include  electron-neutral  collisions  than  to  the  coilison- 
less  calculations. 
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